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Abstract 

In (2+l)-dimensional QED with a Chern-Simons term, we study dynamical break- 
ing of chiral symmetry, using the Dyson-Schwinger equation for the fermions. There is 
a region in parameter space were dynamical chiral symmetry breaking occurs, just as 
in pure QED3 (without Chern-Simons term); outside this region, this chiral symmetry 
breaking solution does not exists. Our results, both numerically and analytically, show 
that the chiral phase transition is a discontinuous first-order transition. 
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1. Introduction 



Dynamical symmetry breaking in (2+l)-dimensional QED has attracted a lot of attention 
over the last decade0)'EI. From a purely theoretical point of view, it is an interesting and 
useful theory to study dynamical mass generation. It also has some applications, not only 
in in elementary particle physics, but also in condensed matter physics, in connection with 
phenomena occurring in planar surfaces. For particle physics, the main interest in three- 
dimensional theories comes from the fact that they are the high-temperature limit of the 
corresponding four- dimensional theory, and thus play a role in studies of the early universe 
and the electroweak phase transition. 

A natural extension of pure QED3 is to add a Chern-Simons (CS) term for the gauge 
field I, which breaks parity explicitly. Recently, it has been shown that in the presence of 
such a CS term there is a first-order chiral phase transition The CS term generates a 
parity-odd mass term for the fermions, but in addition there might be a parity-even mass, 
which breaks chiralityi- ) '0 ) . Both numerical and analytical studies of the truncated Dyson- 
Schwinger (DS) equation for the fermion propagator show that there is a first-order phase 
transition This is quite remarkable and in contrast to the infinite-order phase transition 
one finds in pure QED30)'&. 

In this paper, based on Ref.[|) in which more details are given, we show how such a first- 
order transition arises, and we also compare some different approximation schemes. In the 
next section, we introduce the formalism we are using. In Sec. § we discuss the DS equation 
and the truncation schemes we are using, and how this approach leads to dynamical mass 
generation in pure QED3. Next, we give the results we find in quenched QED3 with CS 
term. In Sec. |5], we include the effects of the vacuum polarization, and finally we give some 
concluding remarks in Sec. [| 



§2. Formalism 

The Lagrangian in Euclidean space is 

C = 4>(i + - m e - rm )i) + \F* V + \i6e pvp A p d u A p + £ gauge fi xing , (1) 

with the dimensionful parameter 9 determining the strength of the CS term. We use four- 
component spinors for the fermions, and the matrix r is defined in such a way that the term 
m ipTip is odd under a parity transformationS)'&. Also the CS term is odd under a parity 
transformation, the other terms in the Lagrangian are parity invariant. 
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The inverse full fermion propagator can be written as 

S-\p) = A e (p) p'+A (p)rp^-B e (p) - B (p)r. 



(2) 



The functions A(p) and B(p) are scalar functions of p 2 , and their bare values are A e — 1, 
A Q = 0, B e = m e , and B Q = m . For convenience, we will introduce another decomposition, 
using the projection operators x± = |(1 ± r )- Defining A± = A e ± A Q and £>± = £> e ± £> G , 
this leads to 



S(p) 



A+(p) ^ + B+(p) _ A_(p) p'+B.jp) 



(3) 



With a four- dimensional representation for the 7-matrices, we can define chirality similar 
as in four- dimensional QED. Without an explicit mass m e for the fermions, the Lagrangian is 
chirally symmetric, but the parity-even mass m e breaks this symmetry. Note that the other 
mass, m , is chirally invariant. Just as in pure QED, the chiral symmetry can be broken 
dynamically due to nonperturbative effects: starting without a mass term for the fermions 
in the Lagrangian, there might be a mass generated by the dynamics. This can be studied 
using the fermion DS equation with the explicit mass m e equal to zero. We also set the bare 
parity-odd mass, m , equal to zero. 

The gauge boson has not only a transverse and longitudinal part, but also a parity-odd 
part, proportional to e^ p g p , due to the explicit CS term in the Lagrangian. In a general 
covariant gauge, the full gauge boson propagator can be written as 



D T (q 2 ) 
D°(q 2 ) - 



q z 

q 2 - n T {q) 



(q2-nT{q)y + {no{q)-e\q\y 

n°(q)-9\q\ 



(4) 
(5) 
(6) 



(q2-nT(q)y + (IlO(q)-6\q\y 

where i7 T and 77° are the transverse and parity-odd part of the vacuum polarization n^ v . 
These parts can be projected out, and this leads to 

d 3 A; 



n 1 (q) = -2Ne 2 



k —2k-q — 3(k- q) 2 /q 2 x 



n°{q) 




(2tt / 

A + (k + q)A + (k) 



q) 2 A 2 + + B 2 + )(k 2 A\ + B 2 + ) 
d 3 A; 

q- (k + q) X 



+ 



A + ,B, 



A_,B_ 



(7) 



(2tt) 3 
A + (k + q)B + (k) 



{{k + q) 2 A 2 + + Bl)(k 2 A 2 + + B 2 



A_,B_ 



(8) 
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using a bare vertex. Note that the longitudinal part of the vacuum polarization is zero 
because of gauge invariance. 

§3. Dynamical Mass Generation 

In order to study dynamical mass generation for the fermions, we use a truncation of the 
DS equation for the fermion propagator. 

3.1. Dyson-Schwinger Equation 

The set of DS equations forms a nonperturbative, exact set of equations for Green's 
functions, and for the fermion propagator it reads 



The problem in solving this equation is the the unknown full vertex function r u (p, k), and 
the vacuum polarization H ^(q) which appears in the gauge boson propagator. In principle, 
we could write down DS equations for these Green's functions as well, but the DS equation 
for the vertex involves a four-point Green's function, and so on. For practical purposes, we 
have to truncate the set of DS equations. 

In this paper, we consider some different truncation schemes. First, we take the most 
simple truncation, namely replacing the full vertex and the full gauge boson by the bare 
ones, see Sec. This is usually referred to as the quenched ladder approximation. We also 
consider a truncation based on the 1/N expansion. 

3.2. 1/N Expansion 

A popular approximation scheme in QED3 with N fermion flavors is the 1/N expansion!). 
The coupling constant e 2 has the dimension of mass, and we can use the large N limit in 
such a way that e 2 j and the product N e 2 remains fixed: Ne 2 = 8a with a fixed. In 
QED3, the most commonly used truncation scheme of the DS equation is based on this 1/N 
expansion. 

In this truncation scheme, the one-loop vacuum polarization has to be taken into account, 
since it is of order one: there are N fermion loops, and each loop is proportional to 1/N. 
The full vertex is replaced by the bare one, because that is the leading order contribution 
in 1/N; all vertex corrections are of down by a factor of 1/N. This truncation scheme gives 
two sets of two coupled integral equations for A + and B + , respectively for A_ and B_ 



S-\ P ) = S \p)-e 2 J- 



d 3 k 
(2^)3 



jfj,S(k)r u (p, k)D^{p - k) . 



(9) 
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V 



A ± (k) ( [D^ q) - » ) *±*JL + a 4^)± B ± (k) D o {q f^ | , (10) 
d 3 k 



q 2 J q 2 2q 

r, , , ' f 1 



(2tt) 3 fc 2 yl|(fc) + l4(fc) 
B ± (k) (2D T (q) + ^jT 2A ± (k)D°(q)^- j . (11) 

where Z) T and are given by Eqs. (5) and (6); a is the gauge parameter. Note that 
the equations for A + and B + are coupled to the ones for A_ and B_ through the vacuum 
polarization only, at least in this truncation scheme. The quenched ladder approximation 
can be recovered by taking a bare gauge boson propagator. 

Of course, we have to make sure that we satisfy at least approximately the Ward- 
Takahashi identity. One of the consequences of the Ward-Takahashi identity is that the 
vertex renormalization and the wave-function renormalization have to be equal. This means 
that if we are using a bare vertex approximation, we have to make sure that A + = A_ = 1, 
at least approximately. To achieve this, we can use the argument that (formally) the wave- 
function renormalization is 1 + 0(1/ N), see Eq. (10) with e 2 = 8a/N, so to leading order in 
1/N we can neglect the wave- function renormalization. Alternatively, we could use a suit- 
able nonlocal (momentum dependent) gauge function a(q), which makes the wave-function 
renormalization equal to one 

3.3. Pure QED3 

In QED3 without CS term, dynamical breaking of chiral symmetry has been studied 
extensively^' %>. It is known that parity is not broken dynamically 0), so the parity-odd 
functions A Q and B Q are zero. Using either a nonlocal gauge such that A(p) = 1, or just 
neglecting the wave-function renormalization, the set of Eqs. (10) and (11) reduces to one 
equation for the dynamical mass function m(p) = B(p)/A(p) 



a f d 3 k m(k) 



a 



m ^ = wJ^wr^m{ ZDT ^ + ^)- (12 » 

With the one-loop vacuum polarization of bare massless fermions, i7 T = — acq, there 
exists a critical number of fermion flavors, N c , above which there is no chiral symmetry 
breaking; for iV < iV c , there is a dynamically generated fermion mass&'B'. Close to the 
critical number, the value of the mass function at the origin behaves like 

m(0) ~ aexp I = I , (13) 



^N c /N-l i 
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corresponding to an infinite-order phase transition. This behavior is obtained both numer- 
ically and analytically, using bifurcation theory. Recently, this critical behavior has been 
confirmed by numerical calculations including the effects of the generated fermion mass on 
the vacuum polarization©. The critical number N c depends on the truncation and gauge: 
in the Landau gauge, N c = 32/tt 2 , whereas in the nonlocal gauge, one finds N c = 128/(37r 2 ). 
Apart from this difference of a factor of |, the results in Landau and nonlocal gauge agree 
with other. 



§4. Results with a Chern-Simons Term in Quenched QED 

One important difference with pure QED, is that with a CS term, there is no trivial solu- 
tion B± = 0; there is always the chirally-symmetric, parity-odd solution B + (p) = —B_(p) = 
B (p), with B e (p) = 0. This solution can also be found using perturbation theory. Dynami- 
cal chiral symmetry breaking can only occur if there is a solution with B + ^ —B-, in such 
a way that B e (p) = ^(B + (p) + B_(p)) is nonzero. Furthermore, in the ultraviolet region, 
the CS term dominates, leading to a behavior of the mass function at large momenta like 
B±(p) ~ ±#/p, whereas the dynamical mass function in pure QED3 falls off like B(p) ~ 1/p 2 - 

4.1. Numerical Results 

In quenched QED3 with a CS term, the equations for A + and B + decouple from the 
equations for v4_ and _B_, as can be seen from Eqs. (10) and (11) with a bare gauge boson. 
By choosing the Landau gauge, a = 0, we find that A = 1 exactly in pure QED, and for small 
values of 9 (small compared to the mass scale e 2 ), A± ~ 1, so we neglect the wave-function 
renormalization, and consider the equations for B + and B- only. 

Numerically, we find the following solutions: 

• B + (p) = —B-(p), the symmetric solution with B e (p) = 0; 

• B + (0) ~ £_(0) > 0, with B e (p) ^ and B (p) = 0(0). 

The last solution exists for small values of 9 only, and corresponds to dynamical breaking of 
chiral symmetry. For large values of 9, we only find the chirally-symmetric solution. Between 
these two regions there is a critical value of 9 for dynamical chiral symmetry breaking. As 
can be seen from Fig. 1, the second solution, B, decreases for increasing 9, and beyond some 
critical value of 9 we cannot find this solution any longer, but at this critical value, B(0) nor 
the chiral condensate go to zero. This clearly indicates that it is a discontinuous, first-order 
phase transition. 
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Fig. 1. B(0) versus 9, both numerically (diamonds: symmetric solution, and plusses: broken so- 
lution) and analytically (dashed: symmetric, and solid line: broken solution). 



4.2. Analytical Results 

Although the numerical results indicate that it is a first-order phase transition, it is 
very difficult to establish such a discontinuous transition numerically. In order to solve 
the equations analytically, we make an expansion in 9, keeping the linear terms only, and 
perform the angular integrals analytically. Next, we use an Ansatz for the solution based on 
the leading ultraviolet behavior of the mass functions B±(p) ~ ±e 2 9/(ir 2 p), but with a finite 
value at the origin. Qualitatively, this Ansatz behaves quite similar to the numerical solution 
we have found in the previous section. This leads to a boundary condition for M = B + (0) 




This cubic equation in M has either one or three 
solutions for given 9, as can be seen from Fig. 2: 
for small (nonzero) values of 9 there are three so- 
lutions, and there is dynamical breaking of chiral 
symmetry, whereas for large values of 9 there is 
only the chirally symmetric solution. From this 
figure it is clear that it is indeed a discontinuous 
transition, at the critical coupling 9 C the "second" 
solution B does not go to zero, nor does the chiral 
condensate. So both our numerical and analyti- 
cal calculations are in good qualitative agreement, 
and show that there is a first-order chiral phase 



l + ln2) + 



7T 



9 M 



0. 



(14) 
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Fig. 2. Boundary condition for M at 
6 = 0.003 and 9 = 0.004, just below 
and above 9 C . 
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transition. Quantitatively there is about a factor of two difference between our numerical 
and analytical results 



6 c /e 2 = 0.00371 (an.) 6 c /e 2 ~ 0.008 (num.) 
B c / e 2 = 0.0352 B c /e 2 ~ 0.06 , 

caused by the approximations we made in solving the equation for B± analytically. 



§5. Influence of the Vacuum Polarization 

Next, we consider the case of N fermion flavors, and include the effects of the vacuum 
polarization. Since we know already that there is a critical number of fermion flavors for 
dynamical chiral symmetry breaking in pure QED, we expect the same to be true in the 
presence of a (small) CS term. We first use the one-loop vacuum polarization with bare 
massless fermions. 

5.1. One-loop Vacuum Polarization 

With massless fermions, there is no parity-odd part of the vacuum polarization, as can 
be seen from Eq. (8), and the transverse part is simply ilr(?) = — ocq, with a = Ne 2 /8. 
Thus the transverse and parity-odd part of the gauge boson become 

W gmq + af + e 2 )' W q 2 ((q + a) 2 + 9 2 ) ' 1 ' 

Again, the equation for A + and B + decouple from the ones for A_ and 5_, and it should 
be noted that if we have a solution (A + , B + ) of the "plus" equations, then the combination 
(A + , —B + ) is automatically a solution of the "minus" equations. 

We used both the Landau gauge and a nonlocal gauged to keep A± ~ 1. In Landau 
gauge, A± will in general not be equal to one, but we nevertheless set it equal to one, in 
order to be consistent with the bare vertex approximation. Most of the results we present 
here, are obtained in the nonlocal gauge; in Landau gauge the results are quite similar, the 
main difference being a scaling factor of | in N, the number of flavors. 

Numerically, we have found results similar to the quenched case, at least below N c {9 = 0), 
the critical number of fermion flavors in pure QED. For fixed 6, we find a critical N c (8), above 
which there is no chiral symmetry breaking, and also for small N, there is a critical 6 C (N) 
for chiral symmetry breaking. Our numerical results indicate a discontinuous (first-order) 
phase transition, as can be seen from Fig. 3(a). 
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Fig. 3. -B(O) (dashed line) and B(0) (solid line) vs. 9 at N = 2 in nonlocal gauge: (a) numerically 
and (b) analytically. For completeness, we also included B e (0) and B o (0) in (a). 



We also solved the equation analytically, in order to show that the phase transition is 
indeed first order. Some standard approximations lead to the general solution 
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B±(p) = M ±2 F l (a, 



} a_ } -;-p 2 /Mi), 



(16) 



with a boundary condition for M± = B±(0) 

16 



±- 



~7T 2 N 



9n 2 N 



3n 2 N 



M± 2 F 1 (a 
M± 2 F 1 (a 



-a 2 /Ml 



+ ,a_,|;-« 2 /M|), 



(17) 



with a± = | ± jiJ (N c /N — 1); the only difference between the Landau and nonlocal gauge 
is the critical number: N c = 32/n 2 in Landau and 128/(37r 2 ) in nonlocal gauge. 

This boundary condition gives us M as function of iV and 8, and is plotted 
in Fig. 3(b). As can be seen from this figure, there 
is a critical 6 C , beyond which the chirally symmet- 
ric solution is the only one. The generated mass 
B does not go to zero at the critical point, cor- 
responding to a first-order phase transition. Also 
in the other direction, increasing iV at fixed 8, we 
find a discontinuous transition, both numerically 
and analytically i* . So both our numerical and an- 
alytical results are in qualitative agreement with 
each other; also both Landau and nonlocal gauge 
are in good qualitative agreement. 

In Fig. 4 we have shown the critical line in the 
(N, 9) plane. 
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are numerical results in nonlocal 
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5.2. Full Vacuum Polarization 



The above results are obtained by using the one-loop vacuum polarization with massless 
fermions. Now one can of course raise the question what is the influence of the generated 
mass on the vacuum polarization. In pure QED3, the generated mass has no effect on the 
critical number of fermion flavors, nor on the behavior of m(0) near the phase transition 
but it does have a crucial effect on the confining properties of this model and on behavior 
of the mass function in the complex momentum plane©. 

Using the vacuum polarization with the dynamical fermion mass, we have both a trans- 
verse and a parity-odd part of the vacuum polarization, and the equations for the "plus" 
and "minus" functions become coupled. For simplicity, we use the Landau gauge, and put 
A± = 1 by hand, so we have to deal with a set of four coupled nonlinear integral equations 
for B + , B_, II T , n° . Our numerical results show that there is dynamical chiral symmetry 
breaking at small values of 9 and N, just as in the previous case. For large values of 9 and N, 
we only find the chirally symmetric solution. In Fig. 5 we have shown our first results. Sofar, 
the numerical results indicate that the phase transition is again discontinuous, but we need 
more accurate results close to the critical parameters before we can establish a first-order 
phase transition; it might be a very sharp second order transition. More detailed calculations 
of the behavior near the critical line should be done, both numerically and analytically. 
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Fig. 5. -B+(0) (diamonds), -B-(O) (plusses), and B + (0) (squares) vs. 8 at (a) N = 0.5 and (b) 
N =1. 



Finally, we should note that an important feature of vacuum polarization is that the 
parity-odd part is linear in the generated fermion mass B. This linearity turns out to 
be important: because of this, the symmetric solution, which is parity-odd, vanishes for 
vanishing 9. This is in agreement with other analyses, showing that there is no dynamical 
breaking of parity in pure QED30), but in contrast to the results obtained with the one-loop 
vacuum polarization of massless fermions, see Fig. 3. 
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6. Conclusion 



In QED3 with N fermion flavors, we have solved the DS equation in the presence of an 
explicit CS term. Without this CS term, there is dynamical chiral symmetry breaking if the 
number of fermion flavors is below some critical number N c . The chiral phase transition at 
N = N c is of infinite order; the order parameter and the generated mass vanish continuously 
at the phase transition. 

In the presence of an explicit CS term for the gauge field, there can be both a chirally 
symmetric mass, and a mass which breaks the chiral symmetry. Below the critical number 
N c , we find that there is chiral symmetry breaking for small values of 9, just as in pure 
QED3. There is a critical line in the (N, 9) plane, beyond which the chirally symmetric 
solution is the only solution. At this critical line, the order parameter and the generated 
mass do not vanish, corresponding to a first-order phase transition. Both numerically and 
analytically we have found this discontinuous transition both in quenched QED3, and using 
the one-loop vacuum polarization. First results including the generated mass in the vacuum 
polarization indicate the same behavior. 

Note that there is no dynamical breaking of parity. In our analysis, this can be seen from 
the fact that the chirally-symmetric, parity-odd solution vanishes for vanishing 9. One has 
to take into account the effect of the generated fermion mass on the vacuum polarization in 
order to get this result. 
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